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individual Dirac eigenvalue distributions and their correlations can be computed directly from the 
effective chiral Lagrangian. 



The XXV International Symposium on Lattice Field Theory 
July 30-4 August 2007 
Regensburg, Germany 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



F K from Dirac eigenvalues at imaginary density 



G. Akemann 



1. Introduction 

The distribution of individual Dirac operator eigenvalues have become a popular tool since 
their calculation ^ in the chiral Random Matrix Theory that is equivalent to the leading-order ex- 
pression for the QCD partition function in the £ -regime The analytical expressions distinguish 
very clearly between different gauge theories and different sectors of topology, as first shown in 
[g|. This has by now been verified by many different groups using different versions of fermions 
with different levels of chiral symmetry on the lattice. It has also become clear how to derive the 
same expression in the £ -regime of chiral perturbations theory (£#PT) [Qj. 

Individual eigenvalue distributions provide perhaps the most efficient tool to extract one of 
the low energy constants (LEC) in #PT, the infinite-volume chiral condensate E. Here we extend 
this analysis to the second LEC in line, the pion decay constant F K , exploiting the fact that a 
nonvanishing chemical potential /I couples to F K to leading order in the £-expansion [Bt]. This 
method has been first suggested for imaginary isospin chemical potential with two different sets 
of Dirac eigenvalues [g]. The advantage over real jl is that the Dirac operator retains its anti- 
hermiticity, allowing for unquenched simulations without encountering any sign-problems, and 
with greatly reduced computer efforts associated with the computation of the lowest eigenvalues. 
The proposal [H was based on the 2-point spectral correlation function computed from £#PT and 
verified their prediction on quenched and unquenched Lattice data. This was generalised in [^] 
where all spectral correlations where computed analytically from the shown equivalence with a 
corresponding chiral Random two-Matrix Theory (chR2MT) with /ii 2- The advantage here is 
that partial quenching is possible, by setting one of the jity to zero. Hence existing configurations 
with jU = can be used to measure F n . This idea was most recently applied to unquenched chiral 
fermions in [ft]. Here we present first results for individual Dirac eigenvalue distributions. 

In section || we introduce the chR2MT and its corresponding £#PT. Section || presents our 
results in a general setting, which is then illustrated pictorially in the simplest case, the quenched 
isospin densities in section f|. Section ^ gives our conclusions and comments on other results. 

2. RMT and #PT with imaginary chemical potential 

We start by defining the chR2MT for imaginary chemical potentials introduced and solved in 



N f 

2?chR2MT ~ / d<Z>cN> exp [-iVTr (<&*<& + ¥*F)] ]\ AsX[9{ji f ) +m f ] . (2. 1) 



The anti-hermitian Dirac matrix £F is given in terms of two complex, rectangular random matrices 
<E> and *P of size N x (N + v) with Gaussian measure: 



Here v corresponds to fixed gauge field topology in the usual way. In the following we restrict our- 
selves to the case of only two different chemical potentials, = ^1,2, withJVi 2 flavours each. 



i 




(2.2) 
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Referring to ref. [Q] for details, we can write down the corresponding eigenvalue representation: 

%hR2MT = / Ft dx i dy i (x i y i ) v+1 (xj+m 2 fl ) ]] (y 2 
J ° i V /1=1 /2=1 

x A w ({A- 2 })A w ({y 2 })det[/ v (2^x^)]^ WLci ^ +C2y ', (2.3) 

where xi and yj are real positive eigenvalues of the matrices <I> + At 1,2*^, respectively. The constants 
c\ t 2 and d depend on fi\^ (see [^). For later convenience we abbreviate the integrand or joint 
probability distribution function (jpdf) by &({x}, {y}). 

If we take the large-Af limit and identify Nxi — ► VLxi = x, Nm.f\ — > VLnifi =rhf\, 2NpL 2 — > 




WF 2 \i\ = jX 2 , and similarly for the second set A2, the partition function eq. (2.3) becomes identical 
to the corresponding £#PT partition function [Jsj] 



J^ pr = y J[/(A^ / )det[[/] v exp 



Tr-^V [U, B] [f/f , 5] + -VLM/Ct/ + t/ 1 ") 



(2.4) 



F w and E have as source terms chemical potential through the charge matrix B =diag(/ii 1^ , /^Im ), 
and the diagonal mass matrix Mf =diag({/n/i }, {nif2}), respectively. For explicit results for these 
partition functions we refer to [ft]. 

3. Results for individual Dirac eigenvalue distributions 

In the following we first define all density correlations, all individual eigenvalues correlations 
(or gap probabilities), and then express the latter in terms of the former. This inversion relation 
is valid for any theory expressed in terms of Dirac eigenvalues, having a jpdf &{{x}, {y}) that is 
symmetric under exchange of all xi and all eigenvalues separately. This applies to our chR2MT 



eq. (2.3), its equivalent £#PT, or a Lattice QCD partition functions in terms of Dirac eigenvalues. 



All density correlation functions are defined by integrating all but k(l) eigenvalues of ^i(^) 



R k>l ( Xu ..., Xk , yi ,..., yi ) = — / n dXi n dyj^({x},{y}) . (3.1) 

J i=k+i j=l+i 



(N-k)\(N-l)\3T. 



The simplest nontrivial example is the probability density R\ t \(x,y) for finding an eigenvalue of 
Qi\ at x and of £^2 at y. If all eigenvalues of one kind are integrated out one finds back the known 
quantities of the one-matrix theory at [i = [§]. Next we define the following gap probabilities 
that the interval [0, s] is occupied by k eigenvalues and [s, °°) by (N — k) eigenvalues of <2)\ , and that 
the interval [0,f] is occupied by / eigenvalues and [?,°°) by (N — I) eigenvalues of 

E kJ (s, t) = p dxi... dx k J dx k+l ...dx N dx x ... dx t J dy M ...dy N 

x ^({x},{y}) , for k,l = 0,1,..., N . (3.2) 

The simplest example is £0,0 (s,?) to find the intervals [0,^] and [0,t] empty of 2>\- and ^-eigenvalues, 
respectively. Similarly we can define the probability to find the k-th ^-eigenvalue at value x k = s, 
and the Z-th f^2 -eigenvalue at value yi = t, to be 

Pkj (s, t) = k ^\ I j Q dx i--- dx k-i j s dx k+\ ■ ■ ■ dxN J^dyi... dyi^x £ dy l+l . . . dy N 

x ^(xi,...,x k -i,x k = s, Xk+ i,...,x N ,yi,...,yi-uyi =t,x t+l ,...,y N ) , (3.3) 
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where the eigenvalues are ordered x\ < ... < x^ and yi < ... < yrf. The simplest example is 
pii(s,t), the distribution of each first eigenvalue. It is easy to see JIo| ] that all the quantities eq. 



( P-3| ) can be obtained from eq. ( p.2| ) by taking two derivatives, 
d 2 

E k> i(s,t) = k\ I ! (pkj (s, t) — Pk+1,1 (s, t) — Pk+\,i (s, t) + Pk+1,1+1 {s, t) ) . (3.4) 



dsd t 



Here, we define p^j = whenever index k or / is zero. Finally we give an inversion formula ex- 
pressing all gap probabilities, and hence all individual eigenvalue distributions in terms of densities: 

N-kN-l rs ,-t 

Ek,i(s,t) = L ■■ -i / dxi...dxk+i / dyi ...dyi +j Rr k+i j +j) (x u ...,x k+i , y u . . . ,y, +j ). 

i=0 7=0 l -J- J0 J0 

(3.5) 

The derivation []!(]] follows closely the pL = case [0]. Since it is known how to generate all higher 
density correlations R^j from resolvents by inserting additional, auxiliary pairs of fermions and 
bosons e.g. into eq. ( |2.4| ), this relation clarifies how to generate individual eigenvalue distributions 
in this setting from field theory. 

4. Examples 

We discuss in detail the simplest example, the probability p\^(s,t). It follows from the gap 
probability using eq. (3.4): 

d 2 

—E $(s,t) = p hl (s,t) . (4.1) 
We expand Eofi(s,t) to include at most 3-point density correlations as an approximation, 

rs f-t rS ft 

E 00 (s,t) = 1- / dxRifl(x)- / dyR ,i(y) + dx dyRn(x,y) 
Jo Jo Jo Jo 

1 f> 1 f s 

+- dyidy 2 R 0}2 {yi,y 2 ) + - dx\dx 2 R 2a {x\,x 2 ) 

2 Jo 2 Jo 

-- [ dx\dx 2 [ dyR 2 ,i(xi,x 2 ,y) - - [ dx [ dyidy 2 Ri, 2 (x,yi,y 2 ) + ...(4.2) 
2. Jo Jo 2 Jo Jo 

The derivatives eliminate all integrals over one-matrix densities that only depend on s or t: 

pi,i(s,t) = Ri,i(s,t) - ( dxR 2 ,i(x,s,t) - ( dyR L2 (s,t,y) + ... . (4.3) 

Jo Jo 

The leading order term is obviously given by the density Ri i (s,t), as can be clearly seen in figs. |] 
and ||. There, we display the quenched density in the case of imaginary isospin chemical potential 
Hi = —\i 2 = In the microscopic large-Af limit Pi,i(£,y) = lim^t»/?i,i(x = x/N,y =y/N) we 
obtain the following result [0, |l] 



P(M) (•*'?) = P(i,0)(*)P(o,i)(?) -xyJtr + (y,x) yX (x,y) - — j/y ( ^2 J e ~^ 

^(x,y) = f l dtte ±2fl2t2 J v (xt)J v (vt) , jr°(jc,y)= f dtt J v (xt)J v {yt) . (4.4) 
Jo Jo 
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Figure 1: Example quenched density R\ \(s,t) (top left) vs individual eigenvalue distribution p\ i(s,t) (top 
right) at v = and 2fi = 0. 159. The lower plots show corresponding 2D cuts at fixed s = 2. 



Here also the well known one-matrix density appears, 



P(l,0)(*) =P(0,i)(*) = ~ [Jl{x) - J v+ i(x)J v ^{x)] = Jf°(x,x) , 



(4.5) 



see fig. |[ Eq. (M) was derived independently for the chR2MT eq. (23) ^ and prior to that for 
£#PT eq. (^4j) using replicas and the Toda-lattice hierarchy It is displayed in figs. [T] and || 
left for topological charge v = and 1, respectively, including 2-dimensional cuts. Because the 
density is the expectation value R\ i (x,y) ~ ( Trd(S>i —x) Tr5(^2 — y) }> A 4 7^ resolves the delta 
function 8(x — y) that we would obtain at pL = 0, times the one-matrix density eq. (-k5) that we 
give for comparison. 

Next we move to individual eigenvalues. A closed determinantal expression for all higher 
density correlation functions in terms of the same building blocks as in eq. ( |4.4| ) was given in ^ : 



pu{{x},{y}) 



k i 

n*<n^ det 



X ( \x h ,x h ) X (x h ,fj 2 ) 



r/v 



■2 J 4/i 2 * v V, 



8/j 2 



(4.6) 



We can insert these formulas into the expansion eq. ( |4.3[ ), after taking the microscopic limit. The 
result truncated at the given order is plotted in figs. [I] and |3] right. The fact that the truncated sum 
is an approximation is seen from the fact that the individual eigenvalue density becomes negative 
(or diverges when adding higher order terms). For the given values in the figs, this happens above 
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s = t m 4, and we have cut the 3D plots at values below —0. 15. Higher order terms in the expansion 



eq. ( |3.5| ) will keep the individual eigenvalue distribution to be zero for larger values of s and t. From 
our experience with the case ji = |Ql we expect that this expansion converges fast. The next to 
leading order used in the figures gives already a reasonably good approximation. 




Figure 3: Same as fig. 1. for v — 1. The exact zero eigenvalues push the density away from the origin. 
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5. Conclusions and outlook 



We have shown how to derive individual eigenvalues distributions p^j for two sets of real Dirac 
operators with different imaginary chemical potentials pi\ 2- Similar expressions have been derived 
for a single Dirac operator with complex eigenvalues at real pL and compared to Lattice data [|TT|]. 
Both types of }i couple to F K and thus allow to fully determine all LECs in the leading order #PT 
Lagrangian. But only imaginary pL with real Dirac eigenvalues allow to date to perform unquenched 
or partially quenched simulations. The equivalence of £#PT to the chR2MT we mentioned here 
for the density and partition function has been derived very recently for all correlators [12]. 

We have given an effective expansion for the distributions by truncating the sum over 
integrated densities, as was illustrated in our examples. The possibility to derive exact expressions 
(which is possible for real /x) is currently under investigation [[K]]. Our hope is that the results 
presented here will become as useful as previously for /I = 0. 



Acknowledgments 

This work was supported by EPSRC grant EP/D03 1613/1 (G.A.) and EU network ENRAGE 
MRTN-CT-2004-005 616. 



References 

[1] S. M. Nishigaki, P. H. Damgaard and T. Wettig, Phys. Rev. D58 (1998) 087704 [hep-th/9803007]; 
P. H. Damgaard and S. M. Nishigaki, Phys. Rev. D63 (2001) 045012 [hep-th/0006 1 1 1 ] . 

[2] E. V. Shuryak and J. J. M. Verbaarschot, Nucl. Phys. A560 (1993) 306 [hep-th/92 12088]. 

[3] R.G. Edwards , Urs M. Heller, Joe Kiskis and Rajamani Narayanan, Phys. Rev. Lett. 82 (1999) 4188. 
[hep-th/9902117]. 

[4] G. Akemann and P. H. Damgaard, Phys. Lett. B583 (2004) 199 [hep-th/031 1 171]. 

[5] D. Toublan and J.J.M. Verbaarschot, Nucl. Phys. B603 (2001) 343 [hep-th/0012144]. 

[6] P. H. Damgaard, U. M. Heller, K. Splittorff and B. Svetitsky, Phys. Rev. D 72 (2005) 091501 

[hep-lat/0508029]; P. H. Damgaard, U. M. Heller, K. Splittorff, B. Svetitsky and D. Toublan, Phys. 
Rev. D 73 (2006) 074023 [hep-lat/0602030]; Phys. Rev. D 73 (2006) 105016 [hep-th/0604054] . 

[7] G. Akemann and T. Wettig, Phys. Rev. Lett. 92 (2004) 102002 Erratum-ibid. 96 (2006) 029902 
[hep-lat/0308003]; J. C. Osborn and T. Wettig, PoS (LAT2005) 200 [hep-lat/05101 15]. 

[8] G. Akemann, P. H. Damgaard, J. C. Osborn and K. Splittorff, Nucl. Phys. B [hep-th/0609059]. 

[9] T. DeGrand, S. Schaefer, arXiv:0708.1731vl (see also in these proceedings). 

[10] G. Akemann and P. H. Damgaard, unpublished, 2007 

[11] G. Akemann, J. Bloch L. Shifrin and T. Wettig, PoS(Lattice2007)224, these proceedings. 
[12] F. Basile, G. Akemann, unpublished, 2007 



7 



